[k,p,err,P]=fixfp('g',100,1.0e-6,100)